Lienard's System and Smale's Problem 



Valery A. Gaiko ^ 

Department of Mathematics 
Belarusian State University of Informatics and Radioelectronics 
L. Beda Str. 6-4, Minsk 220040, Belarus 



Abstract 

In this paper, using geometric properties of the field rotation parameters, we present 
a solution of Smale 's Thirteenth Problem on the maximum number of limit cycles 
for Lienard's polynomial system. We also generalize the obtained result and present 
a solution of Hilbert's Sixteenth Problem on the maximum number of limit cycles 
surrounding a singular point for an arbitrary polynomial system. Besides, we con- 
sider a generalized Lienard's cubic system with three finite singularities, for which 
the developed geometric approach can complete its global qualitative analysis: in 
particular, it easily solves the problem on the maximum number of limit cycles in 
their different distribution. We give also an alternative proof of the main theorem for 
the generalized Lienard's system applying the Wintner-Perko termination principle 
for multiple limit cycles and discuss some other results concerning this system. 
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1 Introduction 

We consider planar dynamical systems 

x = Pnix,y), y = Qrr{x,y), (1.1) 
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where Pn{x,y) and Qn{x,y) are polynomials with real coefficients in the real 
variables x, y, and, ffist of all, a special case of (1.1): classical Lienard's poly- 
nomial system of the form 

x = y, y = -X + niy + fi2y'^ + f^sy^ + ■ ■ ■ + IJ^2ky'^'' + IJ'2k+iy'^''^^- (1-2) 

The main problem of qualitative theory of such systems is Hilbert's Sixteenth 
Problem on the maximum number and relative position of their limit cycles, 
i.e., closed isolated trajectories of (1.1). This problem was formulated as one 
of the fundamental problems for mathematicians of the XX century, however 
it has not been solved even in the simplest (quadratic, cubic, etc.) cases of 
the polynomial systems. In this paper, we suggest a new geometric approach 
to solving the problem in the case of Lienard's system (1.2). In this special 
case, it is considered as Smale's Thirteenth Problem becoming one of the main 
problems for mathematicians of the XXI century [T7|, 

In Section 2 of this paper, applying a canonical system with field rotation 
parameters and using geometric properties of the spirals ffiling the interior and 
exterior domains of limit cycles, we present a solution of Smale 's Thirteenth 
Problem for Lienard's polynomial system (1.2). In Section 3, by means of 
the same geometric approach, we generalize the obtained result and present a 
solution of Hilbert's sixteenth problem on the maximum number of limit cycles 
surrounding a singular point for an arbitrary polynomial system. In Section 4, 
we consider a generalized Lienard's cubic system with three finite singularities, 
for which the developed geometric approach can complete its global qualitative 
analysis: in particular, it easily solves the problem on the maximum number 
of limit cycles in their different distribution. In this section, we give also an 
alternative proof of the main theorem for the generalized Lienard's system 
applying the Wintner-Perko termination principle for multiple limit cycles 
and discuss some other results concerning this system. 



2 Lienard's polynomial system 

System (1.2) and more general Lienard's systems have been studying in nu- 
merous works (see, for example, P, |2], % [7], [I3], [I5], [l8]-[^, [21]). It 
is easy to see that (1.2) has the only finite singularity: an anti-saddle at the 
origin. At infinity, system (1.2) for k ^ 1 has two singular points: a node at 
the "ends" of the j/-axis and a saddle at the "ends" of the x-axis. For studying 
the infinite singularities, the methods applied in pQ for Rayleigh's and van der 
Pol's equations and also Erugin's two-isocline method developed in [TTj can be 
used. Following [11], we will study limit cycle bifurcations of (1.2) by means 
of a canonical system containing only the field rotation parameters of (1.2). 
It is valid the following theorem. 
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Theorem 2.1. Lienard's polynomial system (1.2) with limit cycles can he 
reduced to the canonical form 

x = y = P, y = -x + fiiy + y^ + fi3y^ + ...+y^'' + fi2k+iy^''''^ = Q, (2.1) 
where /xi, /X3, . . . , fi2k+i o^re field rotation parameters 0/ (2.1). 

Proof. Vanish aU odd parameters of (1.2), 

x = y, ?/= -x + /i2?/^ + /i4l/^ + ... + Ai2fey^'', (2.2) 
and consider the corresponding equation 

dy -X + H2y'^ + fJ^^y^ + . . . + iJ^2k y'^'' 



dx y 



Fix,y). (2.3) 



Since F{x, —y) = —F{x,y), the direction field of (2.3) (and the vector field 
of (2.2) as well) is symmetric with respect to the x-axis. It follows that for 
arbitrary values of the parameters /i2, /i4, . . . , fi2k system (2.2) has a center 
at the origin and cannot have a limit cycle surrounding this point. Therefore, 
without loss of generality, all even parameters of system (1.2) can be supposed 
to be equal, for example, to one: fi2 = fJ"i = ■ ■ ■ = fJ'2k = 1 (they could be also 
supposed to be equal to zero). 

To prove that the rest (odd) parameters rotate the vector field of (2.1), let us 
calculate the following determinants: 

Ami = PQ',. - QPU = > 0, 
= PQ',, - QP',s = 1/' > 0, 

Am.... = PQ',..,. - QP',..,. = y'>o. 

By definition of a field rotation parameter [5], for increasing each of the param- 
eters fii, /is, ... , fi2k+i, under the fixed others, the vector field of system (2.1) 
is rotated in positive direction (counterclockwise) in the whole phase plane; 
and, conversely, for decreasing each of these parameters, the vector field of 
(2.1) is rotated in negative direction (clockwise). 

Thus, for studying limit cycle bifurcations of (1.2), it is sufficient to consider 
canonical system (2.1) containing only its odd parameters, /ii, /is, . . . , fi2k+i, 
which rotate the vector field of (2.1). The theorem is proved. 

By means of canonical system (2.1), let us study global limit cycle bifurcations 
of (1.2) and prove the following theorem. 

Theorem 2.2. Lienard's polynomial system (1.2) has at most k limit cycles. 
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Proof. According to Theorem 2.1, for the study of hmit cycle bifurcations of 
system (1.2), it is sufficient to consider canonical system (2.1) containing only 
the field rotation parameters of (1.2): /^i, /i^, . . . , H2k+i- 

Vanish all these parameters: 

x = y, y = -x + y^ + y^ + ... + t/'^ (2.4) 

System (2.4) is symmetric with respect to the and has a center at the 

origin. Let us input successively the field rotation parameters into this system 
beginning with the parameters at the highest degrees of y and alternating 
with their signs. So, begin with the parameter //2fe+i and let, for definiteness, 
Ii2k+i > : 

x^y, y = -x + y^ + y^ + ... + y^'' + fi2k+i (2.5) 

In this case, the vector field of (2.5) is rotated in positive direction (counter- 
clockwise) turning the origin into a nonrough unstable focus. 

Fix ii2k+i and input the parameter ii2k-i < into (2.5): 

x = y, y = -x + y' + y^ + ... + ^2k-i y''"' + 1/'' + ^2k+i (2.6) 

Then the vector field of (2.6) is rotated in opposite direction (clockwise) and 
the focus immediately changes the character of its stability (since its degree 
of nonroughness decreases and the sign of the field rotation parameter at 
the lower degree of y changes) generating a stable limit cycle. Under further 
decreasing /J2k-i, this limit cycle will expand infinitely, not disappearing at 
infinity (because of the parameter ii2k+i at the higher degree of y) . 

Denote the limit cycle by Fi , the domain outside the cycle by Di , the domain 
inside the cycle by D2 and consider logical possibilities of the appearance of 
other (semi-stable) limit cycles from a "trajectory concentration" surrounding 
the origin. It is clear that, under decreasing the parameter /i2fc-i, a semi-stable 
limit cycle cannot appear in the domain D2, since the focus spirals filling this 
domain will untwist and the distance between their coils will increase because 
of the vector field rotation. 

By contradiction, we can also prove that a semi-stable limit cycle cannot 
appear in the domain Di. Suppose it appears in this domain for some values 
of the parameters fjilk+i > A*2fc-i < 0- Return to initial system (2.4) and 
change the inputting order for the field rotation parameters. Input first the 
parameter //2fe-i < 0: 

x^y, y^-x + y'' + y^ + ... + ii2k-iy^''~^ + y'"'- (2.7) 

Fix it under H2k-i = l^2k-i- The vector field of (2.7) is rotated clockwise 
and the origin turns into a nonrough stable focus. Inputting the parameter 
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/i2fc+i > into (2.7), we get again system (2.6), the vector field of wliicli 
is rotated counterclockwise. Under this rotation, a stable limit cycle Fi will 
immediately appear from infinity, more precisely, from a separatrix cycle of 
the Poincare circle form containing infinite singularities of the saddle and node 
types llj. This cycle will contract, the outside spirals winding onto the cycle 
will untwist and the distance between their coils will increase under increasing 
/i2fe+i to the value /i2fc+i- It follows that there are no values of fi2k-i < 
/^2fc+i > O5 which a semi-stable limit cycle could appear in the domain Di. 

This contradiction proves the uniqueness of a limit cycle surrounding the origin 
in system (2.6) for any values of the parameters fi2k-i and /i2fc+i of different 
signs. Obviously, if these parameters have the same sign, system (2.6) has no 
limit cycles surrounding the origin at all. 

Let system (2.6) have the unique limit cycle Fi. Fix the parameters yU2fc+i > 0, 
fJ'2k-i < and input the third parameter, ^2k-3 > 0, into this system: 

x = y, y=-x + y^ + ... + ii^k-z V^^'^ + Z/'''' + • • • + f^2k+i (2.8) 

The vector field of (2.8) is rotated counterclockwise, the focus at the origin 
changes the character of its stability and the second (unstable) limit cycle, 
F2, immediately appears from this point. Under further increasing fi2k-3, the 
limit cycle F2 will join with Fi forming a semi-stable limit cycle, F12, which will 
disappear in a "trajectory concentration" surrounding the origin. Can another 
semi-stable limit cycle appear around the origin in addition to F12? It is clear 
that such a limit cycle cannot appear neither in the domain Di bounded on 
the inside by the cycle Fi nor in the domain D3 bounded by the origin and F2 
because of increasing the distance between the spiral coils filling these domains 
under increasing the parameter fi2k-3- 

To prove impossibility of the appearance of a semi-stable limit cycle in the 
domain D2 bounded by the cycles Fi and F2 (before their joining), suppose the 
contrary, i.e., for some set of values of the parameters, ft2k+i > 0' f^2k-i < 0' 
and /i2fc_3 > 0, such a semi-stable cycle exists. Return to system (2.4) again 
and input first the parameters fi2k-3 > and /i2fc+i > : 

x = y, y = -x + y' + ... + ^,2k-3y"'-' + y"'~' + y'' + fi2k+iy''^'- (2.9) 

Both parameters act in a similar way: they rotate the vector field of (2.9) 
counterclockwise turning the origin into a nonrough unstable focus. 

Fix these parameters under fi2k-3 = /^2fc-3' f^2k+i = fJ'2k+i input the 
parameter fi2k-i < into (2.9) getting again system (2.8). Since, on our as- 
sumption, this system has two limit cycles for fi2k-i > lAk-ii there exists some 
value of the parameter, /i2fc-i {lAk-i < /^2fc-i < 0)) which a semi-stable 
hmit cycle, F12, appears in system (2.8) and then splits into a stable cycle. 
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Fi, and an unstable cycle, under further decreasing /i2fc-i- The formed 
domain D2 bounded by the limit cycles Fi, r2 and filled by the spirals will 
enlarge since, on the properties of a field rotation parameter, the interior un- 
stable limit cycle r2 will contract and the exterior stable limit cycle Fi will 
expand under decreasing ^2k~i- The distance between the spirals of the do- 
main D2 will naturally increase, what will prohibit from the appearance of a 
semi-stable hmit cycle in this domain for /i2fc-i < fJ'll-i- 

Thus, there are no such values of the parameters, fi2k+i > O5 /^2fc-i < O5 
and fi2k-3 > O5 which system (2.8) would have an additional semi-stable 
limit cycle. Obviously, there are no other values of the parameters ^2k+ii 
/i2fc-i, and /i2fc-3 for which system (2.8) would have more than two limit cycles 
surrounding the origin. Therefore, two is the maximum number of limit cycles 
for system (2.8). This result agrees with [23], where it was proved for the first 
time that the maximum number of limit cycles for Lienard's system of the 
form 

i = y, y = -X + niy + nzv^ + iJ-bV^ (2.10) 

was equal to two. 

Suppose that system (2.8) has two limit cycles, Ti and r2 (this is always 
possible if /i2fc+i > -l^2k-i > yU2fc-3 > 0), fix the parameters /i.2fc+i, l^2k-i, 
H2k-3 and consider a more general system than (2.8) (and (2.10)) inputting 
the fourth parameter, fi2k-5 < 0, into (2.8): 

x = y, y = -x + y^ + ... + fi2k-5 Z/''"' + Z/''"' + • • • + f^2k+i (2.11) 

Under decreasing fi2k-5, the vector field of (2.11) will be rotated clockwise and 
the focus at the origin will immediately change the character of its stability 
generating the third (stable) limit cycle, T3. Under further decreasing /i2fe-5, 
Fa will join with F2 forming a semi-stable limit cycle, F23, which will disappear 
in a "trajectory concentration" surrounding the origin; the cycle Fi will expand 
infinitely tending to the Poincare circle at infinity. 

Let system (2.11) have three limit cycles: Fi, F2, F3. Could an additional 
semi-stable limit cycle appear under decreasing ^2k-5, after splitting of which 
system (2.11) would have five limit cycles around the origin? It is clear that 
such a limit cycle cannot appear neither in the domain D2 bounded by the 
cycles Fi and F2 nor in the domain D4 bounded by the origin and F3 because 
of increasing the distance between the spiral coils filling these domains under 
decreasing fi2k-5- Consider two other domains: Di bounded on the inside by 
the cycle Fi and bounded by the cycles F2 and F3. As before, we will prove 
impossibility of the appearance of a semi-stable limit cycle in these domains 
by contradiction. 

Suppose that for some set of values of the parameters fi2k+i > 0^ /^2fc-i < 
/^2fc-3 > 0) /^2fc-5 < 0) such a semi-stable cycle exists. Return to system 



6 



(2.4) again, input first the parameters /i2fc-5 < 0, /i2fc-i < and then the 
parameter ^2k+i > : 

x = y, y = -x+y^+. . .+fi2k-5y^^~^+- ■ ■+f^2k^iy^'''^+y^^+f^2k+iy^''^^ ■ (2.12) 

Fix the parameters /i2fc-5, /^2fc-i under the values fi2k-5y /^2fc-i5 respectively. 
Under increasing /X2fc+i, the node at infinity will change the character of its 
stability, the separatrix behaviour of the infinite saddle will be also changed 
and a stable limit cycle, Fi, will immediately appear from the Poincare circle at 
infinity [IJ. Fix /i2fc+i under the value fi2k+i input the parameter fi2k-3 > 
into (2.12) getting system (2.11). 

Since, on our assumption, (2.11) has three limit cycles for fi2k-3 < /^2fc-3' 
there exists some value of the parameter (0 < /i2fc_3 < fJ^2k-3) ^'^^ which 

a semi-stable limit cycle, F23, appears in this system and then splits into an 
unstable cycle, F2, and a stable cycle, F3, under further increasing fj,2k-3- The 
formed domain bounded by the limit cycles F2, F3 and also the domain 
Di bounded on the inside by the limit cycle Fi will enlarge and the spirals 
filling these domains will untwist excluding a possibility of the appearance of 
a semi-stable limit cycle there. 

All other combinations of the parameters /i2fc+i, fJ'2k-i, /^2fc-3, and fi2k-5 are 
considered in a similar way. It follows that system (2.11) has at most three 
limit cycles. If we continue the procedure of successive inputting the odd 
parameters, fi2k-7, ■ ■ ■ , /^s, /^i, into system (2.4), it is possible first to obtain 
k limit cycles (/i2fc+i > -/i2fc-i > /i2fc-3 > -^2k-5 > At2fc-7 > • • •) and then 
to conclude that canonical system (2.1) (i.e., Lienard's polynomial system 
(1.2) as well) has at most k limit cycles. The theorem is proved. 



3 An arbitrary polynomial system 

Let us consider an arbitrary polynomial system 

x = Pn{x,y,fii,...,fik), y = Qn{x,y,fii,...,fik) (3.1) 

containing k field rotation parameters, /xi, . . . ,fik, and having an anti-saddle 
at the origin. Generalizing the main result of the previous section on the 
maximum number of limit cycles surrounding a singular point in Lienard's 
polynomial system (1.2), we prove the following theorem. 

Theorem 3.1. Polynomial system (3.1) containing k field rotation parameters 
and having a singular point of the center type at the origin for the zero values 
of these parameters can have at most k — 1 limit cycles surrounding the origin. 
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Proof. Vanish all parameters of (3.1) and suppose that the obtained system 

X = Pn{x, y,0,...,0), y = Q„(a;, 0, . . . , 0) (3.2) 

has a singular point of the center type at the origin. Let us input successively 
the field rotation parameters, /ii, . . . , /Xfc, into this system. 

Suppose, for example, that /ii > and that the vector field of the system 

a; = P„(x,?/,/ii,0, . . . ,0), y = Qnix,y,iJ,i,0,...,0) (3.3) 

is rotated counterclockwise turning the origin into a stable focus under in- 
creasing /ii. 

Fix fii and input the parameter /i2 into (3.3) changing it so that the field of 
the system 

X = P„(a;,?/,/ii,/i2,0, . . . ,0), ?/ = (5„(x,?/,/ii,/i2,0, . . . ,0) (3.4) 

would be rotated in opposite direction (clockwise). Let be so for /i2 < 0. 
Then, for some value of this parameter, a limit cycle will appear in system 
(3.4). There are three logical possibilities for such a bifurcation: 1) the limit 
cycle appears from the focus at the origin; 2) it can also appear from some 
separatrix cycle surrounding the origin; 3) the limit cycle appears from a 
so-called "trajectory concentration". In the last case, the limit cycle is semi- 
stable and, under further decreasing /i2, it splits into two limit cycles (stable 
and unstable), one of which then disappears at (or tends to) the origin and the 
other disappears on (or tends to) some separatrix cycle surrounding this point. 
But since the stability character of both a singular point and a separatrix cycle 
is quite easily controlled [11], this logical possibility can be excluded. Let us 
choose one of the two other possibilities: for example, the first one, the so- 
called Andronov-Hopf bifurcation. Suppose that, for some value of /X2, the 
focus at the origin becomes non-rough, changes the character of its stability 
and generates a stable limit cycle, Fi. 

Under further decreasing /i2, three new logical possibilities can arise: 1) the 
limit cycle Fi disappears on some separatrix cycle surrounding the origin; 
2) a separatrix cycle can be formed earlier than Fi disappears on it, then it 
generates one more (unstable) limit cycle, F2, which joins with Fi forming 
a semi-stable limit cycle, F12, disappearing in a "trajectory concentration" 
under further decreasing /i2; 3) in the domain Di outside the cycle Fi or in 
the domain D2 inside Fi, a semi-stable limit cycle appears from a "trajectory 
concentration" and then splits into two limit cycles (logically, the appearance 
of such semi-stable limit cycles can be repeated). 

Let us consider the third case. It is clear that, under decreasing /X2, a semi- 
stable limit cycle cannot appear in the domain D2, since the focus spirals filling 
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this domain will untwist and the distance between their coils will increase 
because of the vector field rotation. By contradiction, we can prove that a 
semi-stable limit cycle cannot appear in the domain Di. Suppose it appears 
in this domain for some values of the parameters /ij > and /Zg < 0. Return 
to initial system (3.2) and change the inputting order for the field rotation 
parameters. Input first the parameter /i2 < : 

i; = P„(a;,y,/i2,0, . . . ,0), y = Q„(a;, /i2, 0, . . . , 0). (3.5) 

Fix it under fi2 = f^2- The vector field of (3.5) is rotated clockwise and the 
origin turns into a unstable focus. Inputting the parameter /ii > into (3.5), 
we get again system (3.4), the vector field of which is rotated counterclockwise. 
Under this rotation, a stable limit cycle, Fi, will appear from some separatrix 
cycle. The limit cycle Fi will contract, the outside spirals winding onto this 
cycle will untwist and the distance between their coils will increase under 
increasing /ii to the value n^. It follows that there are no values of /ig < and 
fil > 0, for which a semi-stable limit cycle could appear in the domain Di. 

The second logical possibility can be excluded by controlling the stability 
character of the separatrix cycle [TT]. Thus, only the first possibility is valid, 
i.e., system (3.4) has at most one limit cycle. 

Let system (3.4) have the unique limit cycle Fi. Fix the parameters /ii > 0, 
/i2 < and input the third parameter, > 0, into this system supposing that 
/is rotates its vector field counterclockwise: 

X = Pn{x,y,fii,^2,fJ'3,0,---,0), y = Qnix,y,fii,fi2,fJ'3,0,---,0)- (3.6) 

Here we can have two basic possibilities: 1) the limit cycle Fi disappears at the 
origin; 2) the second (unstable) limit cycle, F2, appears from the origin and, 
under further increasing the parameter /is, the cycle F2 joins with Fi forming a 
semi-stable limit cycle, F12, which disappears in a "trajectory concentration" 
surrounding the origin. Besides, we can also suggest that: 3) in the domain D2 
bounded by the origin and Fi, a semi-stable limit cycle, F2S, appears from a 
"trajectory concentration" , splits into an unstable cycle, F2, and a stable cycle, 
Fs, and then the cycles Fi, F2 disappear through a semi-stable limit cycle, F12, 
and the cycle Fs disappears through the Andronov-Hopf bifurcation; 4) a semi- 
stable hmit cycle, F34, appears in the domain D2 bounded by the cycles Fi, 
F2 and, for some set of values of the parameters, fil, /i2, /is, system (3.6) has 
at least four limit cycles. 

Let us consider the last, fourth, case. It is clear that a semi-stable limit cycle 
cannot appear neither in the domain Di bounded on the inside by the cycle 
Fi nor in the domain D3 bounded by the origin and F2 because of increasing 
the distance between the spiral coils filling these domains under increasing the 
parameter /is. To prove impossibility of the appearance of a semi-stable limit 
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cycle in the domain D2, snpposc the contrary, i.e., for some set of values of 
the parameters, ^\ > 0, < 0? ^"^^ t^l > 0) such a semi-stable cycle exists. 
Return to system (3.2) again and input first the parameters /i3 > 0, /xi > 0: 

i; = P„(x,y,/xi,/X3,0, ...,0), y = (5n(x, y, /xi, /X3, 0, . . . , 0). (3.7) 

Fix these parameters under /is = /ig, = fil and input the parameter /X2 < 
into (3.7) getting again system (3.6). Since, on our assumption, this system 
has two limit cycles for 112 > A*2) there exists some value of the parameter, 

< < 0), for which a semi-stable limit cycle, ri2, appears in system (3.6) 
and then splits into a stable cycle, Fi, and an unstable cycle, F2, under further 
decreasing /X2- The formed domain D2 bounded by the limit cycles Fi, F2 and 
filled by the spirals will enlarge, since, on the properties of a field rotation 
parameter, the interior unstable limit cycle F2 will contract and the exterior 
stable limit cycle Fi will expand under decreasing /X2- The distance between 
the spirals of the domain D2 will naturally increase, what will prohibit from 
the appearance of a semi-stable limit cycle in this domain for 112 < IJ^- 

Thus, there are no such values of the parameters, /x^ > 0, /i2 < 0, > 0, for 
which system (3.6) would have an additional semi-stable limit cycle. Therefore, 
the fourth case cannot be realized. The third case is considered absolutely 
similarly. It follows from the first two cases that system (3.6) can have at 
most two limit cycles. 

Suppose that system (3.6) has two limit cycles, Fi and F2, fix the parameters 
111 > 0, /i2 < 0, /i3 > and input the fourth parameter, ii^ < 0, into this 
system supposing that /X4 rotates its vector field clockwise: 

= P„(x,y,//i,...,//4,0, ...,0), y = (5„(x, y, //i, . . . , //4, 0, . . . , 0). (3.8) 

The most interesting logical possibility here is that when the third (stable) 
limit cycle, F3, appears from the origin and then, under preservation of the 
cycles Fi and F2, in the domain D3 bounded on the inside by the cycle F3 
and on the outside by the cycle F2, a semi-stable limit cycle, F45, appears and 
then splits into a stable cycle, F4, and an unstable cycle, F5, i. e., when system 
(3.8) for some set of values of the parameters, /xj, fi2^ /ig, /X4, has at least five 
limit cycles. Logically, such a semi-stable limit cycle could also appear in the 
domain Di bounded on the inside by the cycle Fi, since, under decreasing 
/X4, the spirals of the trajectories of (3.8) will twist and the distance between 
their coils will decrease. On the other hand, in the domain D2 bounded on 
the inside by the cycle F2 and on the outside by the cycle Fi and also in the 
domain D4 bounded by the origin and F3, a semi-stable limit cycle cannot 
appear, since, under decreasing fi^, the spirals will untwist and the distance 
between their coils will increase. To prove impossibility of the appearance of a 
semi-stable limit cycle in the domains D3 and Di, suppose the contrary, i.e., 
for some set of values of the parameters, nl > 0, < 0) A'-s > 0) A*4 < 0) 
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such a semi-stable cycle exists. Return to system (3.2) again, input first the 
parameters /i4 < 0, /i2 < and then the parameter /xi > : 



X = P„(x,?/,/ii,/i2,/U4,0, . . . ,0), y = (5„(x,?/,/ii,/i2,/i4,0, . . . ,0). (3.9) 

Fix the parameters /i4, /i2 under the values fil, fi2^ respectively. Under increas- 
ing Hi, a separatrix cycle is formed around the origin generating a stable limit 
cycle, Fi. Fix fii under the value /x^ and input the parameter /is > into (3.9) 
getting system (3.8). 

Since, on our assumption, system (3.8) has three limit cycles for /is < /tg, there 
exists some value of the parameter /t|^ (0 < /t|^ < /ig) for which a semi-stable 
hmit cycle, F23, appears in this system and then splits into an unstable cycle, 
F2, and a stable cycle, Fs, under further increasing /is. The formed domain 
bounded by the limit cycles F2, Fs and also the domain Di bounded on the 
inside by the limit cycle Fi will enlarge and the spirals filling these domains 
will untwist excluding a possibility of the appearance of a semi-stable limit 
cycle there. 

All other combinations of the parameters /ii, fi2, fJ^s, and fi^ are considered in 
a similar way. It follows that system (3.8) has at most three limit cycles. If we 
continue the procedure of successive inputting the field rotation parameters, 
/is, /ig, . . . , /ifc, into system (3.2), it is possible to conclude that system (3.1) 
can have at most k — 1 limit cycles surrounding the origin. The theorem is 
proved. 



4 A generalized Lienard's system 

In [13], we considered a generalized Lienard's cubic system of the form: 

x = y, y = -X + {X - fi)y + (3/2) + fixy - {1/2) + ax'^y . (4.1) 

This system has three finite singularities: a saddle (1, 0) and two antisaddles — 
(0,0) and (2,0). At infinity system (4.1) can have either the only nilpotent 
singular point of fourth order with two closed elliptic and four hj^erbolic 
domains or two singular points: one of them is a hyperbolic saddle and the 
other is a triple nilpotent singular point with two elliptic and two hyperbolic 
domains. We studied global bifurcations of limit and separatrix cycles of (4.1), 
found possible distributions of its limit cycles and carried out a classification 
of its separatrix cycles. We proved also the following theorems. 

Theorem 4.1. The foci of system (4.1) can be at most of second order. 
Theorem 4.2. System (4.1) has at least three limit cycles. 
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Using the results obtained in [T3] and applying the approach developed in this 
paper, we can easily prove a much stronger theorem. 

Theorem 4.3. System (4.1) has at most three limit cycles with the follow- 
ing their distributions: ((1,1), 1), ((1,2),0), ((2, 1),0), ((1,0), 2), ((0,1), 2), 
where the first two numbers denote the numbers of limit cycles surrounding 
each of two anti-saddles and the third one denotes the number of limit cycles 
surrounding simultaneously all three finite singularities. 

Theorem 4.3 agrees, for example, with the earlier results by Iliev and 
Perko [in], but it does not agree with a quite recent result by Dumortier 
and Li [6] published in the same journal. The authors of both papers use very 
similar methods: small perturbations of a Hamiltonian system. In [16], the 
zeros of the Melnikov functions are studied and, in particular, it is proved 
that at most two limit cycles can bifurcate from either the interior or exte- 
rior period annulus of the Hamiltonian under small parameter perturbations 
giving a generalized Lienard system. In [B], zeros of the Abelian integrals are 
studied and it is "proved" that at most four limit cycles can bifurcate from 
the exterior period annulus. Thus, Dumortier and Li "obtain" a configuration 
of four big limit cycles surrounding three finite singularities together with the 
fifth small limit cycle which surrounds one of the anti-saddles. 

The result by Dumortier and Li [6J also does not agree with the Wintner- 
Perko termination principle for multiple limit cycles [11], [23]. Applying the 
method as developed in [3], [8]-[13], we can show that system (4.1) cannot 
have neither a multiplicity-three limit cycle nor more than three limit cycles 
in any configuration. That will be another proof of Theorem 4.3 (the same 
approach can be applied to proving Theorems 2.2 and 3.1 as well). But first let 
us formulate the Wintner-Perko termination principle [23] for the polynomial 
system 

i = /(x,/x), (4.2^) 
where a? G R^; G R"; / G R^ (/ is a polynomial vector function). 

Theorem 4.4 (Wintner— Perko termination principle). Any one-para- 
meter family of multiplicity-m limit cycles of relatively prime polynomial sys- 
tem (4.2^) can be extended in a unique way to a maximal one-parameter 
family of multiplicity-m limit cycles of (4.2^1^) which is either open or cyclic. 

If it is open, then it terminates either as the parameter or the limit cycles be- 
come unbounded; or, the family terminates either at a singular point of (4.2^), 
which is typically a fine focus of multiplicity m, or on a {compound) separatrix 
cycle of {4:.2fj,), which is also typically of multiplicity m. 
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The proof of this principle for general polynomial system {4:.2fx) with a vector 
parameter /x G R" parallels the proof of the planar termination principle for 

the system 

x = P{x,y,X), y = Q{x,y,X) (4.2^) 

with a single parameter A G R (see [11], [23]), since there is no loss of generality 
in assuming that system (4.2^) is parameterized by a single parameter A; i. e., 
we can assume that there exists an analytic mapping /x(A) of R into R" such 
that (4.2^) can be written as (4.2 ^(a)) or even (4.2a) and then we can repeat 
everything, what had been done for system (4.2a) in [23]. In particular, if A is 
a field rotation parameter of (4.2a), it is valid the following Perko's theorem 
on monotonic families of limit cycles. 

Theorem 4.5. If Lq is a nonsingular multiple limit cycle of (4.2o), then Lq 
belongs to a one-parameter family of limit cycles of {4.2 x); furthermore: 

1) if the multiplicity of Lq is odd, then the family either expands or contracts 
monotonically as A increases through Aq; 

2) if the multiplicity of Lq is even, then Lq befurcates into a stable and an 
unstable limit cycle as A varies from Aq in one sense and Lq disappears as A 
varies from \q in the opposite sense; i. e., there is a fold bifurcation at Aq. 

Proof of Theorem 4.3. The proof is carried out by contradiction. Suppose 
that system (4.1) with three field rotation parameters. A, /i, and a, has three 
limit cycles around, for example, the origin (the case when limit cycles sur- 
round another focus is considered in a similar way). Then we get into some 
domain in the space of these parameters which is bounded by two fold bifur- 
cation surfaces forming a cusp bifurcation surface of multiplicity-three limit 
cycles. 

The corresponding maximal one-parameter family of multiplicity-three limit 
cycles cannot be cyclic, otherwise there will be at least one point corresponding 
to the limit cycle of multiplicity four (or even higher) in the parameter space. 
Extending the bifurcation curve of multiplicity-four limit cycles through this 
point and parameterizing the corresponding maximal one-parameter family 
of multiplicity-four limit cycles by a field-rotation parameter, according to 
Theorem 4.5, we will obtain a monotonic curve which, by the Wintner-Perko 
termination principle (Theorem 4.4), terminates either at the origin or on some 
separatrix cycle surrounding the origin. Since we know absolutely precisely at 
least the cyclicity of the singular point (Theorem 4.1) which is equal to two, 
we have got a contradiction with the termination principle stating that the 
multiplicity of limit cycles cannot be higher than the multiplicity (cyclicity) 
of the singular point in which they terminate. 
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If the maximal one-parameter family of multiplicity-three limit cycles is not 
cyclic, on the same principle (Theorem 4.4), this again contradicts to Theo- 
rem 4.1 not admitting the multiplicity of limit cycles higher than two. More- 
over, it also follows from the termination principle that neither the ordinary 
separatrix loop nor the eight-loop cannot have the multiplicity (cyclicity) 
higher than two (in that way, it can be proved that the cyclicity of three 
other separatrix cycles [H] is at most two). Therefore, according to the same 
principle, there are no more than two limit cycles in the exterior domain sur- 
rounding all three finite singularities of (4.1). Thus, system (4.1) cannot have 
neither a multiplicity-three limit cycle nor more than three limit cycles in any 
configuration. The theorem is proved. 

So, we have found two approaches to solving Smale's Thirteenth and Hubert's 
Sixteenth Problems. Both these approaches are based on the application of field 
rotation parameters which determine limit cycle bifurcations of polynomial 
systems. 
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